Coherent population trapping resonances with linearly polarized light 
for all-optical miniature atomic clocks 
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We present a joint theoretical and experimental characterization of the coherent population trap- 
ping (CPT) resonance excited on the Di line of 87 Rb atoms by bichromatic linearly polarized laser 
light. We observe high-contrast transmission resonances (up to ~ 25%), which makes this excitation 
scheme promising for miniature all-optical atomic clock applications. We also demonstrate cancella- 
tion of the first-order light shift by proper choice of the frequencies and relative intensities of the two 
laser field components. Our theoretical predictions are in good agreement with the experimental 
results. 

PACS numbers: 42.72.-g, 42.50.Gy, 32.70.Jz 



I. INTRODUCTION 



There is great interest in the development of minia- 
ture (chip-scale) atomic clocks, magnetometers and other 
metrology tools. One of the promising approaches uti- 
lizes coherent population trapping (CPT) resonances: 
by coherent interaction with two resonant optical fields 
in a A configuration, atoms can be prepared in a 
non-interacting quantum superposition of two hypcrfinc 
ground states (a "dark" state), which results in enhanced 
transmission [lj. The dark state is very sensitive to de- 
tuning of the frequency difference of the two optical fields 
from the hyperfine splitting of the atomic states, which 
allows using the CPT effect for precision sensing. In re- 
cent years atomic frequency standards based on CPT 
resonances with fractional frequency stability of 10~ 12 
or better have been demonstrated P, d, Q • Further im- 
provements in stability require realization of CPT res- 
onances with optimized parameters, such as larger am- 
plitude and contrast and smaller linewidth. In addition, 
any sensitivities to environmental changes must be min- 
imized. In particular, light shifts (i.e., spectra shifts in 
the transmission peak with total light power) can limit 
the frequency stability of CPT-based devices. 

Miniature atomic clocks typically use a single diode 
laser (usually a vertical cavity surface-emitting laser, or 
VCSEL) phase-modulated at an appropriate frequency 
to obtain two light fields, required for CPT resonances. 
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FIG. 1: (Color online) Light-induced two-photon transitions 
in a bichromatic field from F g 
variants of optical polarization. 



As a result, practical operational schemes are restricted 
to light fields of the same polarization. In the case of cir- 
cularly polarized light fields, for example, this restriction 
leads to significant reduction of CPT resonance ampli- 
tude, since large fraction of the atomic population be- 
comes trapped in the "end" Zeeman states with max- 
imum value of the angular momentum (mp = F or 
nip = —F), and not in the desired magnetic field - insen- 
sitive states with mj? = [4. To mitigate this problem, 
a variety of CPT schemes TBI, f| 0, [I| as well as related 
coherent resonances d, OH El have been proposed and 
studied by different groups. 

In this paper we present a CPT interaction scheme that 
takes advantage of the unique level combination of alkali 
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metal atoms with nuclear spin 7=3/2 to realize high- 
contrast magneto-insensitive CPT resonances with two 
optical fields of the same linear polarization {lin\\lin con- 
figuration). This scheme, initially proposed in Ref. [l2l ]. 
is possible when two ground states with total angular 
momentum F g = l,2 are coupled only through an excited 
state with total angular momentum F e =l. In a vapor 
cell this situation is realized only for the Di line of 87 Rb 
when the exited-state hyperfine levels are spectrally re- 
solved. Here, we present an experimental and theoretical 
study of this lin\\lin CPT scheme, and demonstrate the 
cancellation of the first-order light shift by proper choice 
of the frequencies and relative intensities of the two op- 
tical fields, which can reduce the sensitivity of resonance 
parameters to laser intensity fluctuations. 

II. THEORETICAL ANALYSIS 

Observation of CPT resonances requires two optical 
fields Ei,2et~ I (" 1 ' 2 *~' £l ' 2 ^l tuned to two-photon resonance 
with the hyperfine transition of an appropriate atom: 
loi — u>2 = ^hfs- In 87 Rb, for example, the hyperfine 
splitting is Ahfs — 6.835 GHz. For miniature atomic 
clock applications, both fields are produced by a single 
phase-modulated laser (typically a VCSEL) and have the 
same polarization. Traditional CPT clocks use circularly 
polarized laser output [l| . In this case the magnetic field- 
insensitive CPT resonance is formed between mp = 
Zeeman sublevels [(0)-(0) clock transition]. For exam- 
ple, for tr_|_-polarized laser fields E^' 2 ' , the corresponding 
CPT dark state is: 

\dark circ ) = Af x (1) 

r v {1) E (1) 1 

\\F g = l,m = 0) - -f±-^ \F g = 2,m = 0) , 

( 1 2) 

where Vq : are the transition matrix elements for \F g = 
1,2, m = 0) — » \F e ,m — 1), and W is a normalization 
coefficient. The disadvantage of such a scheme is the 
existence of the noninteracting trap state \trap C i rc ) = 
\F g = 2, to = 2): since it is decoupled from both circu- 
larly polarized optical fields, a large portion of the atomic 
population is optically pumped in the trap state. This 
effectively reduces the density of interacting atoms, re- 
sulting in lower CPT resonance contrast. 

Application of linearly (or in general elliptically) polar- 
ized optical fields can eliminate such trap states. How- 
ever, such polarizations simultaneously deteriorate the 
ground-state coherence formed between two the magnetic 
field-insensitive sublevels to = in the presence of a lon- 
gitudinal magnetic field. In this case, two opposite circu- 
larly polarized components of the optical fields simulta- 
neously form two A-systems between the to = ground 
states, through both \F e ,m — ±1) excited states. The 
matrix elements for \F g = 2, to = 0) — > \F e ,m = ±1} 

(2) (2) 

are equal and opposite Vq, = — Vq_ , but those for the 



transitions \F g — 1,to = 0) — > \F e ,m = ±1) are the 

same VqV — Vq_ . As a result, the dark states for each A 
system are different, as can be seen from Eq. [TJ In par- 
ticular, for linearly polarized light, when the intensities 
of the two circular components are the same, the ground- 
state coherence on the (0)-(0) clock transition disappears 
completely, since the dark state for one A system becomes 
the bright (interacting) state for the other. 

87 Rb atoms, with the nuclear spin of 7=3/2, offer 
a unique possibility to form high contrast magneto- 
insensitive CPT dark resonances using a different combi- 
nation of Zeeman sublevels with arbitrary polarized light 
fields even for the condition of reduced or zero coherence 
on the (0)-(0) clock transition. Possible Di optical tran- 
sitions from F g = 1,2 ground state to the F e = 1 excited 
state resulting from the interaction of two elliptically po- 
larized optical fields are shown in Fig. [T] Solid arrows 
show two possible two-photon A resonances that create 
atomic coherence between \F g —l,m——l), \F g —2,m=+l) 
and between \F g =l,m=+l), \F g =2, m=—i) [further in 
the text we refer to these as (— l)-(+l) and (+l)-(— 1) 
resonances] . Both of these A-schemes are excited through 
a common excited state \F e =l,m=0). Most importantly, 
the gyromagnetic g^-factors of the two ground-state hy- 
perfine levels in I = 3/2 alkali metal atoms have equal 
absolute value and opposite signs (neglecting the small 
contribution of nuclear spin) . Thus both ground states in 
each of these A systems shift by an equal amount when a 
weak longitudinal magnetic field is applied. As a result a 
CPT resonance occurs at the unshifted "clock" frequency 
to i — ui2 = ^hfs via two dark states \dark±): 

\dark ± )=Af±x (2) 

r v {1) E {1) 1 

| |F, = l,m = Tl) - \F g = 2, to = ±1) j . 

fl 2) 

Here V±\ ' are the transition matrix elements shown in 
Fig.l, and A/± are normalization constants. In this con- 
figuration there are no "trap" states when the Zeeman 
levels are not degenerate. Hence under CPT conditions 
the atomic population accumulats only in the coherent 
dark states, resulting in higher CPT resonance contrast. 

While the states ([2]) exist for arbitrary elliptical po- 
larizations of the optical fields, a linearly polarized laser 
output is preferable since the influence of nuclear spin 
leads to a small difference in the g-factors for different 
hyperfine components of the ground state. As a result of 
this g-factor difference the frequencies of the (— l)-(+l) 
and (+l)-(— 1) two-photon A-resonances are also differ- 
ent in finite magnetic field (for 87 Rb this difference is 2.8 
kHz/G), but their position is symmetric relative to the 
hyperfine splitting frequency A/j/ s . Thus, for the lin\\lin 
excitation scheme, this small difference in the two-photon 
resonance frequencies leads only to the broadening of the 
resonance, but not to a frequency shift (in linear Zee- 
man effect approximation). At stronger magnetic field 
the CPT resonances from the two A systems become re- 
solvable; the dip between the two CPT peak may be used 
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as a "pseudo-resonance" for atomic clocks as well, al- 
though its characteristics look less promising [l3|, [l4j • 

It is important to note that if the ground states 
F 9 = l,m=±l) and \F g =2,m=±l) are coupled to the ex- 
cited states |F e =2,m=±2), then the dark state given by 
Eq.([2]) are degraded [Hj]. Thus, the key requirement to 
the observation of high-contrast lin\ \ lin CPT resonances 
is good spectral resolution of the excited state F e =l. 
The Dl line of 87 Rb is the most promising candidate 
for this scheme since its excited state hyperfine split- 
ting (812 MHz) exceeds the Doppler broadening at room 
temperature, allowing its use in vapor cells and hence 
in miniature atomic clocks. For all other 7 = 3/2 al- 
kali atoms ( 7 Li, 23 Na, 39 ' 41 K), and also for the D2 line 
of 87 Rb, good spectral resolution of the excited state is 
possible only in the case of laser cooled atoms or in a col- 
limated atomic beam. Even for the 87 Rb Dl line the real- 
ization of a miniature lin\ \ lin CPT atomic clocks will re- 
quire careful optimization of the buffer gas pressure to in- 
crease the diffusion time of atoms in the laser beam and, 
hence, the dark state lifetime without severe degradation 
of the resonance amplitude due to the reduced spectral 
resolution of the pressure-broadened excited states. In 
particular, there is a potential difficulty in using this 
scheme for chip-scale atomic clocks that typically use 
very high buffer gas pressures, which broaden the opti- 
cal resonances beyond the usual Doppler linewidth. This 
difficulty can be circumvented, however, through the ad- 
dition of an anti-relaxation wall coating to a vapor cells 
with lower buffer gas pressure [l6[. 



III. EXPERIMENTAL PROCEDURE 

A schematic of the experimental setup is shown in Fig- 
ure [2] To produce the required resonant bichromatic 
field the output of an external cavity diode laser (ECDL) 
was externally phase-modulated at the 87 Rb hyperfine 
frequency 6.835 GHz by an electro-optical modulator 
(EOM). The laser frequency was stabilized to a satura- 
tion spectroscopy resonance in a vacuum reference cell, 
with the drive field [unmodulated carrier] tuned in the 
vicinity of the 5S\/2F g = 2 <-> hP\i2F e = 1 transition, 
and the probe field [first high-frequency modulation side- 
band] resonant with the 5Si/ 2 F g = 1 <-» 5Pi/ 2 F e = 1 
transition. The two-photon detuning between the probe 
and drive fields was precisely controlled by varying the 
output frequency of the microwave synthesizer providing 
the EOM modulation signal. The intensity ratio I1/I2 
between probe and drive fields [sideband/carrier ratio] 
was adjusted up to 40 % by changing the power of the 
microwave signal sent to the EOM. To achieve higher in- 
tensity ratios we used a thermo-stabilized quartz Fabry- 
Perot etalon with free spectral range of 20 GHz, finesse 
of 30, and maximum transmission of 70%, temperature- 
tuned to partially transmit the strong drive field. The 
beam reflected from the etalon contained the unaffected 
first order sidebands and an attenuated drive field. The 
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FIG. 2: (Color online) Schematic of the experimental setup. 
See text for abbreviations. 



attenuation was controlled by fine-tuning the transmis- 
sion resonance of the etalon via temperature. The re- 
flected beam was collimated to a diameter of 7 mm and 
linearly polarized by a quartz polarizer before entering 
the Rb vapor cell. 

All experimental measurements were recorded using a 
cylindrical Pyrex cell (diameter = 2.5 cm, length = 4 cm) 
containing isotopically enriched 87 Rb and 4 Torr of Ne 
buffer gas to increase the Rb interaction time with the 
laser beam. The cell was mounted inside three-layers of 
magnetic shielding and maintained at a constant temper- 
ature 40° C. To separate the magneto- insensitive "clock" 
resonances (\F g — l,m = —1) «-> \F g = 2, m = +1) 
and \F g = l,m = +1) <-> \F g = 2,m = -1)) from the 
field-dependent resonances, a uniform longitudinal mag- 
netic field of approximately 180 mG was applied using a 
solenoid mounted inside the innermost magnetic shield. 
A photodetector (PD) placed after the cell was used to 
measure the transmitted optical signal. CPT resonance 
parameters such as linewidth, amplitude, etc., were de- 
termined for each set of experimental conditions. 

To continually monitor the CPT resonance frequency 
the EOM modulation frequency was locked to the CPT 
resonance using phase-sensitive detection. A slow fre- 
quency modulation at f m = 400 Hz was superimposed 
on the 6.835 GHz microwave synthesizer output; then 
the photodetector signal was demodulated at f m using a 
lock-in amplifier. The resulting error signal was fed back 
to lock the frequency of the voltage controlled crystal 
oscillator (VCXO) that controlledthe synthesizer. The 
frequency of the locked VCXO was measured by beating 
it with a stable reference derived from a hydrogen maser. 



IV. EXPERIMENTAL RESULTS 

Fig. [Ha) shows typical measured CPT resonances in 
the lin\\lin configuration. The strong peak at zero two- 
photon detuning is from the magneto-insensitive (— 1)- 
(+1) and (+!)-(—!) CPT resonances; the two smaller 



4 



» ( a ) -I 




LU 



(b) 



0.004 



0.002 



0.000 



i 1 r 

-400 -200 200 400 

Two-photon detuning (kHz) 

B=100mG 




~r 

-200 200 

Two-photon detuning (kHz) 

FIG. 3: (Color online) (a) Examples of measured CPT res- 
onances at B = and B = 180 mG. Laser intensity before 
cell is 8.4 mW/cm 2 . (b) Calculated excited-state population 
7r e versus two- photon detuning at B = and B = 100 mG. 
Total resonant intensity h + I2 = 10 mW/cm 2 . 



side peaks correspond do the magnetic field- sensitive res- 
onances formed between the \F g = 2, m = ±2) and \F g = 
l,m = ±0} states. Note that when the Zeeman levels 
are degenerate (at B = 0), coupling to the trap states 
\F g = 2, m = ±2} significantly reduces the amplitude and 
contrast of the CPT resonance (top line). A small mag- 
netic field eliminates coupling to these trap states, result- 
ing in much higher resonance contrast. FigJ3Jb) shows 
the calculated dependence of the total excited state pop- 
ulations 7r e as function of two-photon detuning. These 
calculated spectra can be directly compared to the mea- 
sured CPT resonances, since ir e is proportional to the 
laser field absorption [171 ] ; small discrepancies arise from 
residual magnetic field gradients. 

We modeled each CPT-resonance lineshape as a func- 
tion of the two-photon (Raman) detuning 8r by a gener- 
alized Lorentzian flTl [Lsj : 



I(S R ) = hg + 



B(6 R -6 )j/2 



(7/2) 2 + (S R - S Q r (7/2)2 + {Sr 



<5o) 2 ' 
(3) 

where 7 is the width of the CPT resonance, So is the 
resonance shift, A and B characterize the amplitudes 
of the symmetric and anti-symmetric Lorentzian com- 
ponents, and Ib g is the total transmission away from 
CPT resonance (background transmission). See the Ap- 
pendix for details. A. B, 7, and So are functions of 



the optical field amplitudes and one-photon detuning 
Sl] their values were extracted by fitting experimentally 
recorded or numerically calculated CPT resonances to 
Eq. [3] Fig. 0Ja) shows the calculated variation of A and 
B when the laser frequency is swept across the optical 
transition F g = 2 — > F e = 1 for typical experimental con- 
ditions: one sees that the symmetric Loretzian compo- 
nent A reaches its maximum near the optical resonance; 
whereas the antisymmetric component B is largest for 
Sl ~ ±200 MHz. This resonance asymmetry can lead to 
shifts in measured resonance position if a phase-sensitive 
detection method is used Thus, it is desirable to 

operate a CPT atomic clock under conditions that lead 
to a symmetric CPT resonance. Fig. 0|b) shows that for 
the lin\\lin configuration, CPT resonances become sym- 
metric (B = 0) when the laser frequency is appropriately 
red-detuned from the F e = 1 excited state. We found 
a similar value for the optimal detuning experimentally. 
Note, that this detuning also corresponds to the maxi- 
mum of the the amplitude for the symmetric component 
A, making this operational conditions attractive for clock 
operation. 

Fig. [5](a) shows the measured dependance of the CPT 
resonance width as a function of total laser power, for 
the optimal laser detuning (Sl = —72 MHz) to achieve 
a symmetric resonance. The results of numerical cal- 
culations (solid lines) follow the experimental data quite 
closely; both show the expected linear power broadening. 
Fig. E^b) shows experimental data and numerical calcu- 
lations for the resonance contrast, defined as the ration 
between the resonance amplitude and the off-resonant 
background transmission C = A/Ib g . The contrast in- 
creases rapidly with the laser power, and then saturates 
at C > 25 % for total incident laser power > 1 mW. This 
maximum contrast greatly exceeds that achieved to date 
with traditional CPT clock resonances, C < 4% Q, al- 
beit with relatively large laser power and resonance width 
- 10 kHz. 

The fractional frequency stability of a microwave oscil- 
lator locked to the CPT resonance is proportional to the 
quality figure q — C/7 - the ratio between the resonance 
contrast and width. Notably, q characterizes the achiev- 
able stability for both shot noise-limited and intensity- 
fluctuation-limited atomic clocks [20(. If the frequency 
stability of an atomic clock is limited by photon shot 
noise, the expected Allan deviation, <t(t) is given by [2ll ] : 



1 frje 1 _ 

cr(r) = — . / r 

4 V h g qvo 



1/2 



(4) 



where v§ — A^fs is the atomic reference frequency, e 
is the electron charge, 77 is the photodetector sensitivity 
(measured optical energy per photoelectron), and r is 
the integration time. For the /m||/m CPT resonances in 
our experiment (optimal laser detuning), the maximum 
q ps 3.2 • 10~ 5 /Hz value was measured at w 800 fiW of 
laser power, which corresponds to a resonance contrast 
C ~ 17% and width 7 ~ 5 kHz. This value implies a shot- 
noise-limited frequency stability of cr(r) ~ 2 • 10 _14 r 1 / 2 . 
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FIG. 4: (a) Calculated values of the lineshape parameters A and B [Eq. that describe the symmetric and antisymmetric 
components of a CPT resonance as a function of laser detuning 8l from the F e = 1 level. In these calculations: sideband/carrier 
ratio I1/I2 = 0.5; homogeneous collisionally-broadened optical linewidth 7 = 20 MHz; ground-state relaxation rate P = 500 Hz, 
and the laser beam diameter 6.9 cm. (b) Calculated CPT resonance asymmetry B/A as a function of laser detuning 5l for 
I1/I2 ~ 0.5, 7 = 20 MHz, r = 500 Hz, beam diameter 6.9 cm. CPT resonance becomes symmetric (B = 0) at 5l — —72 MHz. 
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FIG. 5: (a)lin\\lvn CPT resonance linewidth 7 and (b) the resonance contrast C (amplitude/background), for optimal laser 
detuning to achieve a symmetric resonance. Experimental data points are shown as diamonds; solid lines represent numerical 
calculations. 



In practice, however, the frequency stability will also be 
limited by residual intensity noise that depends strongly 
on the linewidth of the laser [22j, |23j • Recently, a table- 
top prototype atomic clock based on lin\ \ tin CPT reso- 
nances with a short-term fractional frequency stability 
2 x 10 _11 t -1 / 2 has been demonstrated [24j . 



V. CANCELLATION OF FIRST-ORDER LIGHT 
SHIFT 

The light shift of a CPT resonance, i.e., changes in the 
maximum transmission caused by changing laser inten- 
sity [25j | , is another important characterization parameter 
for clock applications. Light shifts deteriorate the stabil- 
ity of clocks, since they directly transfer any instabili- 
ties in the laser intensity into variation of the measured 
clock frequency. In the lin\ \ tin configuration the sensi- 
tivity of the CPT resonance position to laser intensity 
may be strongly reduced by choosing a proper intensity 
ratio between two optical fields. For a resonant A system 
the main contribution to the CPT resonance light shift 



comes from the coupling to non-resonant hyperfine lev- 
els, and in particular to the upper hyperfine level F e =2, 
since S hfs = 812 MHz< A hfs = 6.835 GHz. To mo- 
tivate first-order light shift cancellation we consider the 
two limiting cases when one of the optical fields is much 
stronger than the other. If \Ei/E-2\^>l the dark reso- 
nance position is mostly affected by the light shift of the 
Zeeman sublevel |F g = l,m=±l) caused by the strong E\ 
field. This Zeeman shift is negative and proportional 
to -IV^E^/Shfs, and results in a positive shift of the 
CPT resonance. In the opposite case \E±/E2\'€i 1, the 
main contribution is due to shift of the Zeeman sublevel 
|F ff =2,m=±l) oc -\V^E 2 \ 2 /S hfs , which in this case 
leads to a negative shift of the CPT resonance. Since in 
general the frequency of the CPT peak will be determined 
by light shifts created by both fields, it is reasonable to 
expect that by changing the sideband/carrier ratio we 
can find a point of light shift cancellation. Our calcula- 
tions show that the compensation condition depends only 
on the proper choice of the ratio of the two fields and not 
on the absolute value of either optical component. For 
the optimal one-photon laser detuning 5l — —72 MHz, 
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FIG. 6: Measured CPT resonance frequency shift as a func- 
tion of laser power for different sideband/carrier ratios (shown 
for each trace). Experimental conditions were otherwise the 
same as for Fig. [5] Ratio of 0.73 corresponds to the light 
shift compensation condition (i.e., CPT resonance frequency 
minimally affected by laser power variation). 
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FIG. 7: Calculated optimal sideband/carrier ratio, for which 
light shift of CPT resonance position is minimally dependent 
on laser power, as a function of laser detuning. 



the calculations predict 7i//2=0.79 for light shift com- 
pensation, which is in a good agreement with the ex- 
perimentally determined I1/I2 = 0.73 (see Fig(5]). The 
compensation ratio depends slightly on laser detuning, 
as is shown in FigJT] 



VI. CONCLUSION 

In this paper we presented theoretical and experimen- 
tal characterization of magnetic field-insensitive CPT res- 
onances induced in 87 Rb vapor by linearly polarized light. 
This excitation scheme allowed us to greatly increase the 
contrast of magneto-insensitive CPT resonances. The 
dependance of the resonance amplitude, linewidth, con- 
trast, quality figure and frequency on total laser inten- 
sity was investigated. We also demonstrated light-shift 
compensation by proper choice of the frequencies and in- 
tensity ratio of the resonant optical fields. Our results 
suggest that Un\ \ tin CPT resonances may be good candi- 



dates for miniature atomic clocks. 
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VIII. APPENDIX: DENSITY MATRIX 
CALCULATIONS 

Our numerical calculations employ a model based on 
the standard density matrix approach (see e.g., [l7[). Un- 
der assumptions of low saturation and total collisional de- 
polarization of the excited state, the ground-state density 



matrix a gg obeys the following equations: 



n,. 



(•5) 



where H g is the unity matrix in the ground-state man- 
ifold, n g = 2(21 + 1) is the total number of sub-states 
in the ground state, and ir e is the total population of 
the excited state. The first source term in Eq.©, in- 
versely proportional to the excited-state radiative life- 
time r e , describes isotropic repopulation of the ground- 
state sublevels due to spontaneous decay of the excited 
states. The second source term, proportional to the 
ground-state relaxation rate T, describes the entrance 
of unpolarized atoms into the laser beam due to diffu- 
sion and collisions. Hereafter we use the following short- 
hand notations: |e) = \F e , m e ) with m e = —F e , . . . ,F e , 
and \i,m) = m) with m = —Fi,...,Fi (i — 1,2). 
The non-Hermitian ground-state effective Hamiltonian 
has the form: 

H e « = ~^^(|l,m)(l,m|-|2,m)(2,m|) (6) 



+ H B + (R) v -i-U g . 

Here 5r = ^1 — ^2 — fs is the two-photon (Raman) de- 
tuning; Hb is the magnetic Hamiltonian, describing the 
linear Zeeman splitting; the angle brackets (. . .) v denote 
averaging over atomic velocities with a Maxwell distribu- 
tion function. The two-photon excitation matrix for an 
atom moving with a velocity v can be written as 
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™ (^L — kv — UJ e ) + ^7/2 , (Ol — + W,' — UJi — LO e ) + 27/2 

m 



where d is the optical transition dipolc moment opera- 
tor, 5l is the laser (one-photon) detuning, 7 is the optical 
linewidth, and Huj e is the energy of the excited-state HF 
levels. Note that Sl and uj e are measured from a com- 
mon zero level (hereafter from the HF level with mini- 
mal momentum F e = I — J e = 1) . In equations © and 
([7]) the Doppler shift kv (k = [ki + &2]/2) and strong 
velocity-changing collisions are taken into account. The 
matrix R contains resonant (first summation) as well as 
off-resonant (second summation) contributions to light 
shifts and optical pumping rates (Hermitian and anti- 
Hermitian parts, respectively). The non-diagonal (i 5^ j) 
elements of the resonant term induce Raman coherence 
between the hyperfine levels of the ground state, which 
is responsible for the CPT dark resonance. 

Here we consider the steady-state regime, setting 
{d/dt)d gg = in ([SJ. As a spectroscopic signal, we 
take the total excited-state population 7r e , which is pro- 
portional to the total light absorption in optically thin 
media or to the total fluorescence. The following pro- 
cedure is used to find 7r e . From l[5]). the ground-state 



density matrix a gg is expressed in terms of ir e , and then 
7r e is calculated from the normalization condition (0) . In 
the general case of arbitrary magnetic field this algebraic 
problem can be solved numerically. 

At intermediate magnetic fields, when the magnetically 
sensitive resonances are well separated from the magneti- 
cally insensitive resonances, and simultaneously when the 
two magnetically insensitive resonances (— l)-(+l) and 
(+l)-(— 1) are well overlapped (i.e. when their width 
is much greater than their Zeeman splitting), the reso- 
nance lincshapc is described with a good accuracy by the 
generalized Lorentzian [l7], [H| : 

Ci(7/2) 2 C 2 (6 R - £0)7/2 

^ (7/2) 2 + (S R - So) 2 (7/2) 2 + (Sr - So) 2 ' 

(8) 

The coefficients in ([SJ (Cj, 7, and So) can be analyti- 
cally expressed in terms of the parameters of the effec- 
tive ground-state Hamiltonian © in a way similar to the 
method used in [TtJ ■ Alternatively, they can be extracted 
from numerical results by a fitting procedure. 
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We present a joint theoretical and experimental characterization of the coherent population trap- 
ping (CPT) resonance excited on the Di line of 87 Rb atoms by bichromatic linearly polarized laser 
light. We observe high-contrast transmission resonances (up to » 25%), which makes this excitation 
scheme promising for miniature all-optical atomic clock applications. We also demonstrate cancella- 
tion of the first-order light shift by proper choice of the frequencies and relative intensities of the two 
laser field components. Our theoretical predictions are in good agreement with the experimental 
results. 



PACS numbers: 42.72.-g, 42.50.Gy, 32.70.Jz 

I. INTRODUCTION 

There is great interest in the development of minia- 
ture (chip-scale) atomic clocks, magnetometers and other 
metrology tools. One of the promising approaches uti- 
lizes coherent population trapping (CPT) resonances: 
by coherent interaction with two resonant optical fields 
in a A configuration, atoms can be prepared in a 
non-interacting quantum superposition of two hypcrfinc 
ground states (a "dark" state), which results in enhanced 
transmission [1]. The dark state is very sensitive to de- 
tuning of the frequency difference of the two optical fields 
from the hyperfine splitting of the atomic states, which 
allows using the CPT effect for precision sensing. In re- 
cent years atomic frequency standards based on CPT res- 
onances with fractional frequency stability of 10 -12 or 
better have been demonstrated [1-3]. Further improve- 
ments in stability require realization of CPT resonances 
with optimized parameters, such as larger amplitude and 
contrast and smaller linewidth. In addition, any sensi- 
tivities to environmental changes must be minimized. In 
particular, light shifts (i.e., spectra shifts in the transmis- 
sion peak with total light power) can limit the frequency 
stability of CPT-based devices. 

Miniature atomic clocks typically use a single diode 
laser (usually a vertical cavity surface-emitting laser, or 
VCSEL) phase-modulated at an appropriate frequency 
to obtain two light fields, required for CPT resonances. 
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FIG. 1: (Color online) Light-induced two-photon transitions 
in a bichromatic field from F g = 1, 2 to F e = 1 for different 
variants of optical polarization. 



As a result, practical operational schemes are restricted 
to light fields of the same polarization. In the case of cir- 
cularly polarized light fields, for example, this restriction 
leads to significant reduction of CPT resonance ampli- 
tude, since large fraction of the atomic population be- 
comes trapped in the "end" Zeeman states with max- 
imum value of the angular momentum (mp = F or 
nip = —F), and not in the desired magnetic field - insen- 
sitive states with mp = [4]. To mitigate this problem, 
a variety of CPT schemes [5-8] as well as related coher- 
ent resonances [9-11] have been proposed and studied by 
different groups. 

In this paper we present a CPT interaction scheme that 
takes advantage of the unique level combination of alkali 



2 



metal atoms with nuclear spin 7=3/2 to realize high- 
contrast magneto-insensitive CPT resonances with two 
optical fields of the same linear polarization (lin\\lin con- 
figuration). This scheme, initially proposed in Ref. [12], 
is possible when two ground states with total angular 
momentum F g =l,2 are coupled only through an excited 
state with total angular momentum F e =l. In a vapor 
cell this situation is realized only for the Di line of 87 Rb 
when the exited-state hyperfine levels are spectrally re- 
solved. Here, we present an experimental and theoretical 
study of this lin\\lin CPT scheme, and demonstrate the 
cancellation of the first-order light shift by proper choice 
of the frequencies and relative intensities of the two op- 
tical fields, which can reduce the sensitivity of resonance 
parameters to laser intensity fluctuations. 



II. THEORETICAL ANALYSIS 

Observation of CPT resonances requires two optical 
fields Ei ) 2e[ -l ('" 1 ' 2t-fel ' 2Z )l tuned to two-photon resonance 
with the hyperfine transition of an appropriate atom: 
ciJi — t^2 = A/i/s- In 87 Rb, for example, the hyperfine 
splitting is Ahfs — 6.835 GHz. For miniature atomic 
clock applications, both fields are produced by a single 
phase-modulated laser (typically a VCSEL) and have the 
same polarization. Traditional CPT clocks use circularly 
polarized laser output [1]. In this case the magnetic field- 
insensitive CPT resonance is formed between rap = 
Zeeman sublevels [(0)-(0) clock transition]. For exam- 
ple, for cr + -polarized laser fields E+' 2 \ the corresponding 
CPT dark state is: 



\dark C i rc ) = N x 
\\F g = l,m = 0) 



(1) 



T/(l) p(l) 

T/(2)p(2) 



F g = 2, m = 0) 



(1 2) 

where V : arc the transition matrix elements for \F g — 
1,2, m = 0) — » \F e ,m = 1), and J\f is a normalization 
coefficient. The disadvantage of such a scheme is the 
existence of the noninteracting trap state \trap C i rc ) = 
\F g = 2, to = 2): since it is decoupled from both circu- 
larly polarized optical fields, a large portion of the atomic 
population is optically pumped in the trap state. This 
effectively reduces the density of interacting atoms, re- 
sulting in lower CPT resonance contrast. 

Application of linearly (or in general clliptically) polar- 
ized optical fields can eliminate such trap states. How- 
ever, such polarizations simultaneously deteriorate the 
ground-state coherence formed between two the magnetic 
field-insensitive sublevels m = in the presence of a lon- 
gitudinal magnetic field. In this case, two opposite circu- 
larly polarized components of the optical fields simulta- 
neously form two A-systems between the to = ground 
states, through both \F e ,m — ±1) excited states. The 
matrix elements for \F g = 2, m = 0) — > \F e ,m = ±1) 



transitions \F g = l,m = 0) — > \F e ,m — ±1) arc the 
same Vq+ = . As a result, the dark states for each A 
system are different, as can be seen from Eq. 1. In par- 
ticular, for linearly polarized light, when the intensities 
of the two circular components are the same, the ground- 
state coherence on the (0)-(0) clock transition disappears 
completely, since the dark state for one A system becomes 
the bright (interacting) state for the other. 

87 Rb atoms, with the nuclear spin of 7=3/2, offer 
a unique possibility to form high contrast magneto- 
insensitive CPT dark resonances using a different combi- 
nation of Zeeman sublevels with arbitrary polarized light 
fields even for the condition of reduced or zero coherence 
on the (0)-(0) clock transition. Possible Di optical tran- 
sitions from F g — 1,2 ground state to the 7^ = 1 excited 
state resulting from the interaction of two elliptically po- 
larized optical fields are shown in Fig. 1. Solid arrows 
show two possible two-photon A resonances that create 
atomic coherence between |F 9 =l,m=— 1), \F g =2,m=+l) 
and between \F g = l, m=+l), \F g =2,m=— 1) [further in 
the text we refer to these as (— l)-(+l) and (+l)-(— 1) 
resonances] . Both of these A-schemes are excited through 
a common excited state \F e —l,m=0). Most importantly, 
the gyromagnetic (/^-factors of the two ground-state hy- 
perfine levels in 7 = 3/2 alkali metal atoms have equal 
absolute value and opposite signs (neglecting the small 
contribution of nuclear spin) . Thus both ground states in 
each of these A systems shift by an equal amount when a 
weak longitudinal magnetic field is applied. As a result a 
CPT resonance occurs at the unshifted "clock" frequency 
u)\ — u>2 = Ahfs via two dark states \dark±): 



\dark±) = A/± x 



(2) 



\F„ 



1,to = =pl) 



,(2) P (2) 



VX>E : 



2, to = ±1) 



Ti 



(1 2) 

Here V± { > are the transition matrix elements shown in 
Fig.l, and J\f± are normalization constants. In this con- 
figuration there are no "trap" states when the Zeeman 
levels are not degenerate. Hence under CPT conditions 
the atomic population accumulats only in the coherent 
dark states, resulting in higher CPT resonance contrast. 

While the states (2) exist for arbitrary elliptical po- 
larizations of the optical fields, a linearly polarized laser 
output is preferable since the influence of nuclear spin 
leads to a small difference in the g-factors for different 
hyperfine components of the ground state. As a result of 
this ^-factor difference the frequencies of the (— l)-(+l) 
and (+l)-(— 1) two-photon A-resonances are also differ- 
ent in finite magnetic field (for 87 Rb this difference is 2.8 
kHz/G), but their position is symmetric relative to the 
hyperfine splitting frequency A/,/ s . Thus, for the lin\\lin 
excitation scheme, this small difference in the two-photon 
resonance frequencies leads only to the broadening of the 
resonance, but not to a frequency shift (in linear Zee- 
man effect approximation). At stronger magnetic field 
the CPT resonances from the two A systems become re- 
solvable; the dip between the two CPT peak may be used 
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Vn_ , but those for the 
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as a "pseudo-resonance" for atomic clocks as well, al- 
though its characteristics look less promising [13, 14]. 

It is important to note that if the ground states 
\F g = \ .m=±.l) and \F g =2,m=±l) are coupled to the ex- 
cited states \F e =2, m=±2), then the dark state given by 
Eq.(2) are degraded [15]. Thus, the key requirement to 
the observation of high-contrast lin\ \ lin CPT resonances 
is good spectral resolution of the excited state F e =l. 
The Dl line of 87 Rb is the most promising candidate 
for this scheme since its excited state hyperfine split- 
ting (812 MHz) exceeds the Doppler broadening at room 
temperature, allowing its use in vapor cells and hence 
in miniature atomic clocks. For all other 7 = 3/2 al- 
kali atoms ( 7 Li, 23 Na, 39 ' 41 K), and also for the D2 line 
of 87 Rb, good spectral resolution of the excited state is 
possible only in the case of laser cooled atoms or in a col- 
limated atomic beam. Even for the 87 Rb Dl line the real- 
ization of a miniature lin\ \ lin CPT atomic clocks will re- 
quire careful optimization of the buffer gas pressure to in- 
crease the diffusion time of atoms in the laser beam and, 
hence, the dark state lifetime without severe degradation 
of the resonance amplitude due to the reduced spectral 
resolution of the pressure-broadened excited states. In 
particular, there is a potential difficulty in using this 
scheme for chip-scale atomic clocks that typically use 
very high buffer gas pressures, which broaden the opti- 
cal resonances beyond the usual Doppler linewidth. This 
difficulty can be circumvented, however, through the ad- 
dition of an anti-relaxation wall coating to a vapor cells 
with lower buffer gas pressure [16]. 



III. EXPERIMENTAL PROCEDURE 

A schematic of the experimental setup is shown in Fig- 
ure 2. To produce the required resonant bichromatic 
field the output of an external cavity diode laser (ECDL) 
was externally phase-modulated at the 87 Rb hyperfine 
frequency 6.835 GHz by an electro-optical modulator 
(EOM). The laser frequency was stabilized to a satura- 
tion spectroscopy resonance in a vacuum reference cell, 
with the drive field [unmodulated carrier] tuned in the 
vicinity of the 5S\/2F g = 2 <-> hP\i2F e = 1 transition, 
and the probe field [first high-frequency modulation side- 
band] resonant with the 5Si/ 2 F g = 1 <-» 5Pi/ 2 F e = 1 
transition. The two-photon detuning between the probe 
and drive fields was precisely controlled by varying the 
output frequency of the microwave synthesizer providing 
the EOM modulation signal. The intensity ratio I1/I2 
between probe and drive fields [sideband/carrier ratio] 
was adjusted up to 40 % by changing the power of the 
microwave signal sent to the EOM. To achieve higher in- 
tensity ratios we used a thermo-stabilized quartz Fabry- 
Perot etalon with free spectral range of 20 GHz, finesse 
of 30, and maximum transmission of 70%, temperature- 
tuned to partially transmit the strong drive field. The 
beam reflected from the etalon contained the unaffected 
first order sidebands and an attenuated drive field. The 
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FIG. 2: (Color online) Schematic of the experimental setup. 
See text for abbreviations. 



attenuation was controlled by fine-tuning the transmis- 
sion resonance of the etalon via temperature. The re- 
flected beam was collimated to a diameter of 7 mm and 
linearly polarized by a quartz polarizer before entering 
the Rb vapor cell. 

All experimental measurements were recorded using a 
cylindrical Pyrex cell (diameter = 2.5 cm, length = 4 cm) 
containing isotopically enriched 87 Rb and 4 Torr of Ne 
buffer gas to increase the Rb interaction time with the 
laser beam. The cell was mounted inside three-layers of 
magnetic shielding and maintained at a constant temper- 
ature 40° C. To separate the magneto- insensitive "clock" 
resonances (F = l,m = —1 <-> F = 2, m = +1 
and F — 1, m = +1 «-> F = 2, m = —1) from the 
field-dependent resonances, a uniform longitudinal mag- 
netic field of approximately 180 mG was applied using a 
solenoid mounted inside the innermost magnetic shield. 
A photodetector (PD) placed after the cell was used to 
measure the transmitted optical signal. CPT resonance 
parameters such as linewidth, amplitude, etc., were de- 
termined for each set of experimental conditions. 

To continually monitor the CPT resonance frequency 
the EOM modulation frequency was locked to the CPT 
resonance using phase-sensitive detection. A slow fre- 
quency modulation at f m = 400 Hz was superimposed 
on the 6.835 GHz microwave synthesizer output; then 
the photodetector signal was demodulated at f m using a 
lock-in amplifier. The resulting error signal was fed back 
to lock the frequency of the voltage controlled crystal 
oscillator (VCXO) that controlledthe synthesizer. The 
frequency of the locked VCXO was measured by beating 
it with a stable reference derived from a hydrogen maser. 



IV. EXPERIMENTAL RESULTS 

Fig. 3(a) shows typical measured CPT resonances in 
the lin\\lin configuration. The strong peak at zero two- 
photon detuning is from the magneto-insensitive (— 1)- 
(+1) and (+!)-(—!) CPT resonances; the two smaller 
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FIG. 3: (Color online) (a) Examples of measured CPT res- 
onances at B = and B = 180 mG. Laser intensity before 
cell is 8.4 mW/cm 2 . (b) Calculated excited-state population 
7r e versus two- photon detuning at B = and B = 100 mG. 
Total resonant intensity h + I2 = 10 mW/cm 2 . 



side peaks correspond do the magnetic field- sensitive res- 
onances formed between the \F g = 2, m = ±2) and \F g = 
l,m = ±0} states. Note that when the Zeeman levels 
are degenerate (at B = 0), coupling to the trap states 
\F g = 2, m = ±2} significantly reduces the amplitude and 
contrast of the CPT resonance (top line). A small mag- 
netic field eliminates coupling to these trap states, result- 
ing in much higher resonance contrast. Fig. 3(b) shows 
the calculated dependence of the total excited state pop- 
ulations 7r e as function of two-photon detuning. These 
calculated spectra can be directly compared to the mea- 
sured CPT resonances, since 7r e is proportional to the 
laser field absorption [17]; small discrepancies arise from 
residual magnetic field gradients. 

We modeled each CPT-resonance lineshape as a func- 
tion of the two-photon (Raman) detuning 8r by a gener- 
alized Lorentzian [17-19]: 



I(S R ) = hg + 



B(6 R -6 )j/2 



(7/2) 2 + (S R - S Q r (7/2)2 + {Sr 



<5o) 2 ' 
(3) 

where 7 is the width of the CPT resonance, So is the 
resonance shift, A and B characterize the amplitudes 
of the symmetric and anti-symmetric Lorentzian com- 
ponents, and Ib g is the total transmission away from 
CPT resonance (background transmission). See the Ap- 
pendix for details. A. B, 7, and So are functions of 



the optical field amplitudes and one-photon detuning 
Sl] their values were extracted by fitting experimentally 
recorded or numerically calculated CPT resonances to 
Eq. 3. Fig. 4(a) shows the calculated variation of A and 
B when the laser frequency is swept across the optical 
transition F g = 2 — > F e = 1 for typical experimental con- 
ditions: one sees that the symmetric Loretzian compo- 
nent A reaches its maximum near the optical resonance; 
whereas the antisymmetric component B is largest for 
Sl ~ ±200 MHz. This resonance asymmetry can lead to 
shifts in measured resonance position if a phase-sensitive 
detection method is used [20]. Thus, it is desirable to 
operate a CPT atomic clock under conditions that lead 
to a symmetric CPT resonance. Fig. 4(b) shows that for 
the lin\\lin configuration, CPT resonances become sym- 
metric (B = 0) when the laser frequency is appropriately 
red-detuned from the F e = 1 excited state. We found 
a similar value for the optimal detuning experimentally. 
Note, that this detuning also corresponds to the maxi- 
mum of the the amplitude for the symmetric component 
A, making this operational conditions attractive for clock 
operation. 

Fig. 5(a) shows the measured dependance of the CPT 
resonance width as a function of total laser power, for 
the optimal laser detuning (Sl = —72 MHz) to achieve 
a symmetric resonance. The results of numerical cal- 
culations (solid lines) follow the experimental data quite 
closely; both show the expected linear power broadening. 
Fig. 5(b) shows experimental data and numerical calcu- 
lations for the resonance contrast, defined as the ration 
between the resonance amplitude and the off-resonant 
background transmission C = A/Ib g . The contrast in- 
creases rapidly with the laser power, and then saturates 
at C > 25 % for total incident laser power > 1 mW. This 
maximum contrast greatly exceeds that achieved to date 
with traditional CPT clock resonances, C < 4% [1], al- 
beit with relatively large laser power and resonance width 
- 10 kHz. 

The fractional frequency stability of a microwave oscil- 
lator locked to the CPT resonance is proportional to the 
quality figure q — C/7 - the ratio between the resonance 
contrast and width. Notably, q characterizes the achiev- 
able stability for both shot noise-limited and intensity- 
fluctuation-limited atomic clocks [21]. If the frequency 
stability of an atomic clock is limited by photon shot 
noise, the expected Allan deviation, <j(t) is given by [22]: 



1 frje 1 _ 

cr(r) = — . / r 

4 V h g qvo 



1/2 



(4) 



where v§ — A^fs is the atomic reference frequency, e 
is the electron charge, 77 is the photodetector sensitivity 
(measured optical energy per photoelectron), and r is 
the integration time. For the lin\\lin CPT resonances in 
our experiment (optimal laser detuning), the maximum 
q ps 3.2 • 10~ 5 /Hz value was measured at w 800 fiW of 
laser power, which corresponds to a resonance contrast 
C ~ 17% and width 7 ~ 5 kHz. This value implies a shot- 
noise-limited frequency stability of cr(r) ~ 2 • 10 _14 r 1 / 2 . 
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FIG. 4: (a) Calculated values of the lineshape parameters A and B [Eq. (3)] that describe the symmetric and antisymmetric 
components of a CPT resonance as a function of laser detuning 8l from the F e = 1 level. In these calculations: sideband/carrier 
ratio h/h = 0.5; homogeneous collisionally-broadened optical linewidth 7 = 20 MHz; ground-state relaxation rate Y = 500 Hz, 
and the laser beam diameter 6.9 cm. (b) Calculated CPT resonance asymmetry B/A as a function of laser detuning 5l for 
I1/I2 = 0.5, 7 = 20 MHz, r = 500 Hz, beam diameter 6.9 cm. CPT resonance becomes symmetric (B = 0) at Sl = —72 MHz. 
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FIG. 5: (a)lin\\lin CPT resonance linewidth 7 and (b) the resonance contrast C (amplitude/background), for optimal laser 
detuning to achieve a symmetric resonance. Experimental data points are shown as diamonds; solid lines represent numerical 
calculations. 



In practice, however, the frequency stability will also be 
limited by residual intensity noise that depends strongly 
on the linewidth of the laser [23, 24]. Recently, a table- 
top prototype atomic clock based on lin\\lin CPT reso- 
nances with a short-term fractional frequency stability 
2 x 10 _11 t -1 / 2 has been demonstrated [25]. 



V. CANCELLATION OF FIRST-ORDER LIGHT 
SHIFT 

The light shift of a CPT resonance, i.e., changes in the 
maximum transmission caused by changing laser inten- 
sity [26] , is another important characterization parameter 
for clock applications. Light shifts deteriorate the stabil- 
ity of clocks, since they directly transfer any instabili- 
ties in the laser intensity into variation of the measured 
clock frequency In the lin\ \ lin configuration the sensi- 
tivity of the CPT resonance position to laser intensity 
may be strongly reduced by choosing a proper intensity 
ratio between two optical fields. For a resonant A system 
the main contribution to the CPT resonance light shift 



comes from the coupling to non-resonant hyperfine lev- 
els, and in particular to the upper hyperfine level F e =2, 
since S hfs = 812 MHz< A hfs = 6.835 GHz. To mo- 
tivate first-order light shift cancellation we consider the 
two limiting cases when one of the optical fields is much 
stronger than the other. If \Ei/E 2 \^>l the dark reso- 
nance position is mostly affected by the light shift of the 
Zeeman sublevel |.F ff =l,m=±l) caused by the strong E\ 
field. This Zeeman shift is negative and proportional 
to —\V^Ei\ 2 /Shf s , and results in a positive shift of the 
CPT resonance. In the opposite case \Ei/E 2 \^_ 1, the 
main contribution is due to shift of the Zeeman sublevel 
|Fg=2,m=±l) oc -{V^ E 2 \ 2 /Shfs, which in this case 
leads to a negative shift of the CPT resonance. Since in 
general the frequency of the CPT peak will be determined 
by light shifts created by both fields, it is reasonable to 
expect that by changing the sideband/carrier ratio we 
can find a point of light shift cancellation. Our calcula- 
tions show that the compensation condition depends only 
on the proper choice of the ratio of the two fields and not 
on the absolute value of either optical component. For 
the optimal one-photon laser detuning Sl = —72 MHz, 
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FIG. 6: Measured CPT resonance frequency shift as a func- 
tion of laser power for different sideband/carrier ratios (shown 
for each trace). Experimental conditions were otherwise the 
same as for Fig. 5. Ratio of 0.73 corresponds to the light 
shift compensation condition (i.e., CPT resonance frequency 
minimally affected by laser power variation). 
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FIG. 7: Calculated optimal sideband/carrier ratio, for which 
light shift of CPT resonance position is minimally dependent 
on laser power, as a function of laser detuning. 



the calculations predict I\f 72=0.79 for light shift com- 
pensation, which is in a good agreement with the ex- 
perimentally determined J1//2 = 0.73 (see Fig. 6). The 
compensation ratio depends slightly on laser detuning, 
as is shown in Fig. 7. 



VI. CONCLUSION 

In this paper we presented theoretical and experimen- 
tal characterization of magnetic field-insensitive CPT res- 
onances induced in 87 Rb vapor by linearly polarized light. 
This excitation scheme allowed us to greatly increase the 
contrast of magneto-insensitive CPT resonances. The 
dependance of the resonance amplitude, lincwidth, con- 
trast, quality figure and frequency on total laser inten- 
sity was investigated. We also demonstrated light-shift 
compensation by proper choice of the frequencies and in- 
tensity ratio of the resonant optical fields. Our results 



suggest that lin\ \ lin CPT resonances may be good candi- 
dates for miniature atomic clocks. 
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VIII. APPENDIX: DENSITY MATRIX 
CALCULATIONS 

Our numerical calculations employ a model based on 
the standard density matrix approach (see e.g., [17]). Un- 
der assumptions of low saturation and total collisional de- 
polarization of the excited state, the ground-state density 
matrix a gg obeys the following equations: 



Tr{a 93 } = 1, 



H e fia gg - d gg H cji 



TTe \ n„ 

— + r 

T e / n a 



(5) 



where H g is the unity matrix in the ground-state man- 
ifold, n g = 2(21 + 1) is the total number of sub-states 
in the ground state, and n e is the total population of 
the excited state. The first source term in Eq.(5), in- 
versely proportional to the excited-state radiative life- 
time r e , describes isotropic repopulation of the ground- 
state sublevels due to spontaneous decay of the excited 
states. The second source term, proportional to the 
ground-state relaxation rate T, describes the entrance 
of unpolarized atoms into the laser beam due to diffu- 
sion and collisions. Hereafter we use the following short- 
hand notations: |e) = \F e , m e ) with m e = —F e , . . . , F e , 
and \i,m) = m) with m = —F l ,...,F i (i = 1,2). 
The non-Hermitian ground-state effective Hamiltonian 
has the form: 

H e s = -y^(|l,m)(l,m|-|2,m)(2,m|) (6) 

m 

+ H B + (R) v -i^U g . 

Here 5r ; = uji — u 2 — A h f s is the two-photon (Raman) de- 
tuning; Hb is the magnetic Hamiltonian, describing the 
linear Zeeman splitting; the angle brackets (. . .)„ denote 
averaging over atomic velocities with a Maxwell distribu- 
tion function. The two-photon excitation matrix for an 
atom moving with a velocity v can be written as 
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R- V \i ^ ^' m l( d ' E ') l e )( e l( d - E j)l^ m ') / 7 - m n | V It m) ( '' m|(d 



where d is the optical transition dipole moment opera- 
tor, Sl is the laser (one-photon) detuning, 7 is the optical 
linewidth, and hu> e is the energy of the excited-state HF 
levels. Note that Sl and uo e are measured from a com- 
mon zero level (hereafter from the HF level with mini- 
mal momentum F e = I — J e = 1). In equations (6) and 
(7) the Doppler shift kv (k — [fci + &2]/2) and strong 
velocity-changing collisions are taken into account. The 
matrix R contains resonant (first summation) as well as 
off-resonant (second summation) contributions to light 
shifts and optical pumping rates (Hcrmitian and anti- 
Hermitian parts, respectively). The non-diagonal (i ^ j) 
elements of the resonant term induce Raman coherence 
between the hyperfine levels of the ground state, which 
is responsible for the CPT dark resonance. 

Here we consider the steady-state regime, setting 
{d/di)d gg = in (5). As a spectroscopic signal, we 
take the total excited-state population 7r e , which is pro- 
portional to the total light absorption in optically thin 
media or to the total fluorescence. The following pro- 
cedure is used to find 7r e . From (5), the ground-state 



density matrix d gg is expressed in terms of 7r e , and then 
7r e is calculated from the normalization condition (5). In 
the general case of arbitrary magnetic field this algebraic 
problem can be solved numerically. 

At intermediate magnetic fields, when the magnetically 
sensitive resonances are well separated from the magneti- 
cally insensitive resonances, and simultaneously when the 
two magnetically insensitive resonances (— 1)— (+1) and 
(+l)-(— 1) are well overlapped (i.e. when their width 
is much greater than their Zccman splitting), the reso- 
nance lineshapc is described with a good accuracy by the 
generalized Lorentzian [17, 19]: 

Ci(7/2) 2 C 2 (6 R - £0)7/2 

^ ° (7/2) 2 + (Sr - So) 2 (7/2) 2 + (Sr - So) 2 ' 

(8) 

The coefficients in (8) (Ci, 7, and i5o) can be analyti- 
cally expressed in terms of the parameters of the effec- 
tive ground-state Hamiltonian (6) in a way similar to the 
method used in [17]. Alternatively, they can be extracted 
from numerical results by a fitting procedure. 
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